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Teaching Mathematics with DERIVE & DERIVE in Education 

Ich habe noch einige -wenige -Exemplare von Teaching Mathematics with DERIVE, Tagungsband 
von Krems 1992 vorratig. Sie konnen das Buch bei mir um oS 330.-, bzw. DM 50.- bestellen. 

Wir haben es auch wieder ubernommen, fur unsere deutschsprachigen Freunde den Tagungsband von 
Krems 93 ’’DERIVE in Education” (siehe Book Shelf) zu importieren. Das Buch ist ausgezeichnet 
gelungen und bietet einen hervorragenden Querschnitt iiber die Einsatzmoglichkeiten von DERIVE. 
Leider ist dieses Buch etwas teurer. Wir konnen es Ihnen um oS 450.-, bzw. DM 65.- anbieten. Die 
Preise verstehen sich inklusive Porto und Verpackung. Mit der Auslieferung des zweiten Buches, 
miissen Sie sich aber bis August gedulden. Jetzt kommt die Plymouth Konferenz, dann ein wenig 
Urlaub und dann sind wir wieder im Einsatz. Anruf, FAX oder ein kurzes Schreiben geniigen. 
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LETTER OF THE EDITOR 


Liebes DUG - Mitglied, 

gerade noch rechtzeitig vor den groBen Ferien wird 
der DNL# 14 fertig. Herzlichen Dank fur die vielen 
Zuschriften und Beitrage. Ich bin sehr erfreut, dass 
die DUG in vermehrtem MaBe kleinere Dienst- 
leistungen fur die Mitglieder verrichten konnte. 
Viele von Ihnen werden sich im nachsten Monat 
anlasslich der DERIVE-Konferenz in Plymouth 
zurn ersten Mai oder wieder treffen. Ich freue mich 
schon darauf, eine groBe Anzahl unserer Mitglieder 
kennen zu lemen, bzw. wieder zu sehen. Das 
Programm verspricht auch einen spannenden 
Tagungsverlauf. 

Mit Dr. S. Biryukov aus Moskau hat die DUG nun 
auch ein Mitglied in Russland. Dr. Biryukov hat 
einen sensationellen Beitrag geschickt: einen 

Zeichensatz fur DERIVE, mit dem sich von 
DERIVE aus Grafiken bezeichnen lassen. Dazu gibt 
es eine Utility zum Generieren eigener 
Zeichensatze. Zwei Screenshots konnen Sie auf 
dieser Seite bewundern. Weiteres folgt im DNL#15. 

Die Auswahl der Artikel fur den DNL erfolgt 
primar nach technischen Gesichtspunkten (Umfang, 
Fort-setzung, ...) und eventuell nach 
zusammengehorigen Themengruppen. Ich bitte die 
Autoren um Geduld. nichts geht verloren. Ich denke 
auch daran, im nachsten Jahr den DNL 
umfangreicher zu gestalten. Die Materialfulle lasst 
dies zu. 

Im nachsten DNL mochte ich einige DERIVE- 
unterstiitzte Schularbeiten osterreichischer Kollegen 
vorstellen: Aufgaben zu dynamischen Systemen und 
Schularbeiten im Teamwork! ! ! 

Mit den besten Wiinschen fur einen schonen 
Sommer 

Ihr 
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Dear DUG Member, 

just in time before summer vacations will start 
DNL# 14 has been finished. Many thanks for lots of 
letters and contributions. It is fine that we can 
support our members with little services more and 
more. Many of you will meet next month again or 
the first time at the DERIVE-Conference in 
Plymouth, UK. I am looking forward to meeting a 
lot of our members and seeing many of you again. 
The conference programme promises some thrilling 
days, too. 

With Dr. S. Biryukov from Moscow the DUG has 
now its first Russian member. He has delivered a 
really sensational contribution: A font set for labe- 
ling and describing DERIVE plots - in the 
DERIVE-environment - together with a tool to 
create own DERIVE fonts. You can find two 
DERIVE-screens on this page ... and wait for 
DNL#15, please! 

Selection of contributions for the DNL is primarily 
a technical question (size, sequel,) and sometimes I 
am looking for a possible connection between the 
articles. I beg the authors for patience, no 
contribution is lost. I have the plan to give the DNL 
more content next year. The plenty of material will 
allow it. 

With my best wishes for a fine summer 


Sincerely yours 



Scale x:2 


Derive ZD- plot 






P 2 


EDITORIAL 
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The DERIVE-NEWSLETTER is the 
Bulletin of the DERIVE User Group. It is 
published at least four times a year with a 
contents of 30 pages minimum. The goals 
of the D-N-L are to enable the exchange of 
experiences made with DERIVE as well as 
to create a group to discuss the possi- 
bilities of new methodical and didactical 
manners in teaching mathematics. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/Fax: +43-(0)2275/8207 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to 
reinforce the international touch of the 
D-N-L. It must be said, though, that non- 
English articles will be warmly welcomed 
nonetheless. Your contributions will be 
edited but not assessed. By submitting 
articles the author gives his consent for 
reprinting it in D-N-L. The more contri- 
butions you will send, the more lively and 
richer in contents the DERIVE Newsletter 
will be. 


Preview: (Contributions for the next issues): 

Stability od systems of ODEs, Kozubik, SLK 

Algebraic Operations on Polynomials in DERIVE, Roanes, ESP 

Prime Iterating Number Generators, Wild, UK 

Graphic Integration, Probability Theory, Linear Programming, Bohm, AUT 
DERIVE in Austrian Schools, some examples, Lechner & Eisler, AUT 
Tilgung fremderregter Schwingungen, Klingen, GER, 

Der Fermat-Punkt im Dreieck, Geyer, GER 
Continued Fractions in DERIVE, Cordoba a.o., ESP 
Turtle-Commands in DERIVE, Lechner, AUT 
Dreieck. MTH, Wadsack, AUT 

Newton Method and Ill-Conditioned Problems, Lopes, POR 

2D Plots Labeling, Biryukov, RUS 

Plotting t-periodic functions, Verhoosel, NED 

IMP Logo and Misguided Missiles, Sawada, HAWAII 

Discussion of Curves, Kuenzer & Reichel & Bohm, IT & AUT 

Reichel-Splines versus Bohm-Splines, a Competition, AUT 

3D-Geometry, Reichel, AUT 

Parallel- and Central Projection, Bohm, AUT 

Conic Sections, Fuchs, AUT 


and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag.Josef Bohm 
Herstellung: Selbstverlag 
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K. Herdt, Osnabriick, Germany 

Approximation/Interpolation 

Bei der Interpolationsaufgabe wird zu vorgegebenen Punkten (Xj,yi), i = 0,1,..., k das Interpolations- 
polynom y=a 0 + a^ + ...+ a k x k bestimmt mit P(Xi) = yj. Bekanntlich fuhren die Verfahren Newton-Inter- 
polation, Lagrange-lnterpolation und der direkte Ansatz uber das lineare Gleichungssystem zum 
selben Ergebnis. Zu Zeiten der Tischrechner (und der Handrechnung) wurde das Newton-Verfahren 
der dividierten Differenzen vorgezogen, da unter anderem leichter die Anzahl der zu berucksich- 
tigenden Punkte geandert werden konnte. Mit Einfuhrung der programmierbaren Rechner wurde die 
Lagrange-lnterpolation bedeutsamer, da ihr Algorithmus in hoheren Programmiersprachen 
wesentlich leichter umgesetzt werden kann als der der Newton-Interpolation. Mit den Computer- 
Algebra-Syste-men wird nun der (direkte) Losungsweg uber das lineare Gleichungssystem 
zunehmend interessant, da die auftretende Vandermonde-Koeffizientenmatrix mit den Mitteln der 
CAS leicht aufgestellt werden kann und die Losung linearer Gleichungssysteme ohnehin integriert ist. 
Hinzu kommt, dass hierbei die Interpolation deutlicher als Grenzfall der Approximationsaufgabe 
I(p(Xi) - yi) 2 = Min. in Erscheinung tritt. Mit den Fehlergleichungen 


Va - y = 5 

(V = Vandermonde-Matrix der x r Potenzen, a = Vektor der gesuchten Polynomkoeffizienten und 
y = Vektor der yi) lasst sich die Approximationsaufgabe auch kurz mit 5 T -5 =Min. beschreiben und a 
ergibt sich bekanntlich als Losung des linearen Gleichungssystems 


V T V a = V T y 


(*)• 


Im Falle einer quadratischen Matrix V wird naturlich 5 = 0 und die Fehlergleichungen bilden gerade 
das lin. Gleichungssystem der Interpolationsaufgabe, brauchten also gar nicht uber (*) gelost zu 
werden. Lediglich aus Praktikabilitatsgrunden wird jedoch auch in diesem Falle (*) verwendet, denn 
dann braucht beim Rechengang nicht zwischen Interpolation und Approximation unterschieden zu 
werden. Die schlechtere Kondition von (*) gegenuber V a = y braucht in diesem Falle nicht zu storen, 
da sich ja die Rechengenauigkeit jederzeit erhohen lasst. Eine „Mustersitzung“ habe ich beigefugt. 


From approximation to interpolation using the 
Vandermonde matrix. Using this method you can 
see that the interpolation problem is the extreme 
case of the approximation problem 
I(p(Xi) - yi) 2 = Min. 

With V a - y = 5 we obtain 5 T 5 =Min. and a is the 
solution of the system of linear equations 
V T V a = V T y. 

(V = Vandermonde matrix of the powers of Xj, a = 
vector of the polynomial coefficients, y = vector of 
the yi) 

#1: Datenmatrix der (xi,yi) - The Matrix of 



the data 


f 1 11 


2 2 


#2: d 


3 4 

4 9 


5 15 


L 8 0 J 


#3: Vandermonde_Matrix arid RHS (k = order of the polynomial): 
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11 ■ x 149 

#31: P0L(d , 1) = + 

37 37 

2 

2485- x 23497- x 4297 

#32: P0L(d , 2) = + 

2904 2904 484 

3 2 

7885-x 65765 ■ x 652211 -x 82727 

#33: P0L(d , 3) = + + 

19868 14901 59604 9934 

4 3 2 

3959 -x 23959 -x 476053 -x 44944 -x 19955 

#34: P0L(d , 4) = + + 

44336 22168 133008 8313 11084 

#35: The interpolation polynomial: 

#36: INT_P0L(d) := P0L(d, DIMENSION(d) - 1) 

5 4 3 2 

41- x 23- x 1705- x 1847- x 20107-x 146 

#37: INT_P0L(d) = + + 

2520 56 504 168 1260 21 


#38: much more elegant: instead of functions KOEFF, P0T_X, P0L(d,k): 

k i 

#39: P0LY(d, k) := I ELEMENT (S0LUTI0N(d , k) , i + 1, k + 2)-x 

i=0 

3 2 

7885- x 65765- x 652211-x 82727 

#40: P0LY(d , 3) = + + 

19868 14901 59604 9934 

#41: INT_P0LY(d) := P0LY(d, DIMENSION(d) - 1) 

5 4 3 2 

41- x 23- x 1705- x 1847-x 20107-x 146 

#42: INT_P0LY(d) = + + 

2520 56 504 168 1260 21 


Additional note of the Editor: 

There is a contribution which is related to the topic above in DNL#63: Beyond Polynomial 
Regression. 

H. Scheuermann. Hofheim/Taunus. Germany 

(I try to give a brief summary of H. Scheuermann's letter) 

I . Letter: ... It would be nice if I could copy one DERIVE expression from one Algebra Window to 
another one (to change the expression or to try another way of solving the problem or ...) 

2. Letter (some days after): ... By chance I've found a way to copy expressions using Ctrl+U. 

Open two Algebra Windows and load two files into the windows 
Highlight the expression to be copied in window #1 
AUTHOR and F3 (or F4) - copying the expression into the Author line 
Change the window pressing FI (the expression seems to be lost but ...) 

AUTHOR (in the other window) and press Ctrl+U - the expression appears in the Edit line ready 
for proceeding as you want. 


DNL: Thanks for the hint, and many thanks for the keyboard model. 
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(1) It is a bit amazing that it is so difficult to simplify expressions only containing numbers, + 1 

and powers [it is easy to construct immensely more complicated examples than (1); in fact, (1) is the 
simplest I can figure out]. Since there are some very clever people at Soft Warehouse, maybe 
someone there [or anyone else, for that matter] can tell us a little about simplifying algorithms; what 
can be implemented and what can't. 


(2) In teaching the nature of mathematics, examples of this kind, in combination with DERIVE, could 
serve as a starting point for a discussion on what a proof is in mathematics. Can we prove that x = 1 
by approximating it to more and more significant digits? One could relate this to "constructive 
mathematics". 

Has anyone on this BBS any experience or thoughts on this? 


In general, we have a tendency on this BBS to discuss the user interface of DERIVE, which is fine, 
but I would find it interesting to discuss also other experiences of DERIVE. For instance, Roger 
Folsom's discussion about the chain rule on this BBS was very useful to me. Jerry Glynn asked some 
time ago if we were interested in discussing teaching with DERIVE. I should have answered YES, but 
didn't since I have no such experience, so I would have little to contribute, BUT, I do teach 
mathematics, and I would be very interested to "listen" to such a discussion. 


Message 3082:From JERRY GL YNN to HARALD LANG about #3081 I SIMPLIFY AND DERIVE 

Wonderful!! I will ignore your disclaimer and take your entire statement as a contribution to a 
discussion about teaching math. I will try out your provocative example and then comment. 


Message 3089: From CHRISAR to PUBLIC about QUADRATIC EQUATIONS AND DERIVE 

Recently some classmates and I were assigned solve some simple quadratic equations using Derive 
in our college's math lab. The following equation caused Derive to produce only one of two correct 
values for x; "x A (2/3) - 2x A (1/3) - 35 = 0." In actuality, x = {-125, 343}, but Derive will only provide 343 
as a possible value. Can anyone tell me why this is, and what this means in terms of Derive's 
reliability? Thanks! 


Message 3090: From JERRY GLYNN to CHRISAR about #3089 / QUADRATIC EQUATIONS AND 
DERIVE 

Your question connects to complex and real numbers. If you plot the left side of your equation you'll 
get a graph which only appears when x > 0. Derive (and many other programs) thinks that (-125) A (2/3) 
is a complex number. In this case -25/2+25 sqrt(3)*i/2. The complex i can be entered by holding down 
the ALT key and typing I ... you get an i with a cap on it. In Derive Manage Branch will get you to a 
choice of Principal or Real or Any. Solve your equation with the setting on Principal (the default) and 
then try it with the setting at Real. Also graph the left side of your equation for each of these settings. 
I'll be happy to discuss this further and respond to any questions or comments. Should you question 
the results you get from Derive? Yes! Also question the answers from all other sources (including me). 
Thanks for your good question! 

See how DERIVE 6 is behaving. Josef 

2/3 1/3 

#1: SOLVE(x - 2-x - 35 = 0, x, Real) 

#2: x = 343 

2/3 1/3 

#3: S01VE(x - 2-x - 35 = 0, x) 


#4: 


x = 343 


8 
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#5: 

2/3 

NSOLVE(x 

1/3 

- 2-x -35=0, x, 

Real ) 


#6: 

X = 342.9999923 



#7: 

2/3 

NSOLVE(x 

1/3 

- 2-x -35=0, x, 

-1000, 

0) 

#8: 





#9: 

2/3 

(-125) 

25 25-73-l 

2 2 



#10: 

2/3 

(-125) 

1/3 

2 • (-125) - 35 r - 

105 

15-73-l 


#11: Branch Real 

2/3 1/3 

#12: (-125) - 2- (-125) - 35 i 0 

2/3 1/3 

#13: SOLVE(x - 2-x - 35 = 0, x, Real) 

#14: 


2/3 

#15: (-125) r 25 


Message 3091: From HADUD to HARALD LANG about SIMPLIFYING CUBE ROOTS 

The equation (SQRT(5)+2) A (1/3)-(SQRT(5)-2) A (1/3) = 1 follows immediately from the fact that 

SQRT(5) +\- 2 = (SQRT(5)/2 +\- 1/2) A 3. 

You shouldn't be "amazed" that DERIVE gives up on this expression. The problem is that fractional 
powers are involved here and, as you point out, there is no general algebraic method (other than an 
infinite series) to expand a fractional power of a sum. It is only in highly special cases that 
simplification is possible. Consider the steps DERIVE would have to go through to try and find a 
simplification: 

1) Fractional powers are multivalued functions; which branch is to be used? In our example both 
arguments of the cube roots are positive. Assuming that the principal branch has been selected by 
the user, there is no ambiguity here. 

2) Inspecting the arguments DERIVE would have to recognize that they are so-called "quadratic 
irrationalities", i.e. members of a field (SQRT(i) + j)/k, where i, j, k are integers and i is not a perfect 
square (a field is a group that is closed under addition, multiplication and division). Therefore there 
is a chance (no certainty !) that the arguments can be represented as the third power of a member 
of the same field. In that case the expression would simplify to a number which, again, belongs to 
the same field. 

3) To check out this possibility expand 

((SQRT(5)+u)/v) A 3 = (3 u A 2 + 5)/v A 3 SQRT(5) + (u A 3 + 15 u)/v A 3 
and set up the equations 

3u A 2 + 5 = v A 3 
u A 3 + 1 5u = 2v A 3 

Do these nonlinear equations have a solution with integer u? 

4) Eliminating v by multiplying the first equation by 2 and subtracting the second equation gives 

u A 3 - 6 u A 2 + 15u-10 = 0, 
the only real root of which is u = 1 . Hence v = 2. 

It is clear from this that in the general case of a fractional exponent m/n DERIVE must find an 
integer solution of a n-th degree polynomial in order to simplify the expression. For n > 4 this 
becomes very difficult and time-consuming. Thus it is understandable (and perhaps desirable) that 
it doesn't bother. 
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Message 3098: From HARALD LANG to HADUD about #3091 / SIMPLIFYING CUBE ROOTS 


Very interesting! So what to do is 

(1) We consider the number field { x*sqrt(n) + y } where n is a fixed positive integer which is not a 
perfect square (n = 5 in my example) and x and y are arbitrary rational numbers. 

(2) We check whether each of the numbers sqrt(5) +\- 2 is a perfect cube in this field, which amounts 
to solving a third-degree equation [BTW: when we 'identify coefficients', we actually use that sqrt(n) is 
irrational, a nice classroom problem.] 

I think you have clarified to us in a very nice way what we ask DERIVE to perform when we try to 
simplify my expression. BTW, it seems that we have used the failure of DERIVE to perform a 
simplification to motivate the introduction of the concept 'number field' [now I am a teacher again.] 

Here is another example: 8*SIN(pi/18)*COS(pi/9)*COS(2*pi/9) 

This is equal to 1 (new challenge: prove it!). We can get reasonably convinced by approximating it by 
DERIVE, but DERIVE fails to simplify it to 1. A Bug!! Hello, Soft Warehouse, are you there 
<there...echo...echo> :-) 

No, no; no bug. But what intelligent can be said about this example? Is there again some number field 
around to help us understand what we demand from DERIVE when we try simplifying? 


[To the sysop: If you think the issue if off topic for this BBS, don't hesitate to delete this message, and 
give me a message to shut up.] --Harald 


Message 3099: From JERRY GLYNN to HARALD LANG about #3098 

But why would I tell you to shut up?? I want to hear what comes next! 


This is how DERIVE 6 treats 
the expression now: 


# 1 : 

# 2 : 

# 3 : 

# 4 : 

# 5 : 



Trigonometry Expand 



4-^2-C0S| — 




Trigonometry Collect 



Some more short notes for the DERIVE USER FORUM 

Bob from Homer. NY, USA 

Sorry I'm late with this (membership dues, ed.) -please don’t cancel me - love your DNL issues. 
Thanks, Bob. 

DNL: Don I worry, we didn’t cancel you. And much pleasure with this DNL. 

G. Scheu. Pfinztal, Germany 

DearJosef, do you know the factorization of a 2 + b 2 = (|a| + |b| + 2 V|a b|) (|a| + |b| - 2 V|a b|). DERIVE 
only knows expanding from right hand side to left hand side. 

In one file you can find some functions for calculating the sum of the digits of an number. 

DNL: We will publish this file in one of the next DNLs. 

Dr. S.V.Birvukov. Moscow. Russia 

I have a utility that supports DERIVE 2D Plots labeling and axes drawing, labeling and numbering. 
Dr. Kutzler told me that it will be rather interesting for DERIVE Users and suggested to describe it in a 
short paper for DUG Newsletter . 
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Thomas Weth: A Lexicon of Curves 


(4) 
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Ebene Algebraische und 
T r a n s z e n d e n t e Kurven (4) 

Thomas Weth, Wurzburg, Germany 

Einheitliche Konstruktionen 

In den ersten drei Folgen dieser Reihe wurden die (neben Kegelschnitten) bekanntesten 
und bedeutendsten algebraischen Kurven (zunachst dritter Ordnung) vorgestellt. Im 
Einzelnen handelte es sich um die Kissoide (Folge 1), die Strophoide (Folge 2) und die 
Trisektrix von MacLaurin (Folge 3). Fur jede dieser Kurven wurde neben der Herleitung 
ihrer algebraischen Gleichung auch jeweils eine elementargeometrische punktweise 
Konstruktionsvorschrift angegeben. Obwohl sich diese Erzeugungsweisen sehr stark 
voneinander unterscheiden, habe ich (mehr oder weniger durch Zufall) eine Moglichkeit 
entdeckt, alle drei Kurventypen durch eine einzige Konstruktion herzustellen. Abgesehen 
von der graphischen Darstellung der Kurven war mir DERIVE vor allem bei der Herleitung 
der algebraischen Gleichungen (und der damit verbundenen endgiiltigen Identifizierung) 
der Kurven die entscheidende Hilfe. Allerdings gelang mir trotz DERIVE nicht die 
Klassifizierung alter im Folgenden auftretenden Kurven. Somit stellt diese Folge des 
"Kurvenlexikons" auch eine Aufforderung an die DERIVE- Spezialisten in aller Weit dar, 
in den genannten Fallen die Identifizierung der Kurven zu versuchen (mit der Bitte, mich 
bei Erfolg davon zu unterrichten). 

Uniform Constructions 

In the first three contributions I have introduced the best known algebraic curves (besides the conic 
sections). The constructions of these curves are very different but by chance I found a possibility to 
produce the curves by only one construction. DERIVE was very helpful not only by plotting the 
graphs but also by deriving the algebraic equations (connected with the identification) of the curves. 
Unfortunately I was not able to classify all of the curves which you will find in this series in future 
contributions. So this part of the "Lexicon of Curves" could be a challenge for all DERIVE-specialists 
to try the identification of the curves - and please contact me in case of success. 

Die Konstruktionsvorschrift - The Rule of Construction 
a) Die Inversions am Kreis - The Inversion wrt a Circle 

Unter der Inversion am Kreis versteht man eine ebene Abbildung gemaB folgender Vorschrift: 

Gegeben ist der Kreis k(M,r) und ein Punkt P ^ M. Der Bildpunkt P’ liegt auf der 
Halbgeraden MP und es gilt: MP • MP’ = r 2 . 

Eine elementargeometrische Konstruktion des Bild- 
punktes P’ gelingt (unter anderen) etwa folgendermaBen: 

Liegt P auBerhalb der Kreislinie, so konstruiert man mit 
der Strecke PM als Hypotenuse das rechtwinklige 
Dreieck PMS (vgl. Zeichnung). Nach dem Kathetensatz 
im rechtwinkligen Dreieck gilt dann fur den LotfuBpunkt 

P’: MP • MF = r 2 . 

Liegt P innerhalb der Kreislinie, so konstruiert man das 
Lot durch P auf PM und erhalt S. Dann vervollstandigt 
man das rechtwinklige Dreieck P’MS und erhalt nach 

dem Kathetenstaz wiederum MP • MP' = r 2 . 


Inversion on a circle: P* is the image point of 
P. If P outside of the circle, then construct the 
right triangle PMS over PM and P’ is the 
pedal point of S wrt to PM and we have 
MP • MP' = r 2 . If P inside of the circle then find 
the perpendicular line through P wrt to PM 
giving intersection point S with the circle and 
accomplish the right triangle MSP’ with the 
right angle in S. Then we have again 
MP-MP = r 2 . 
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a) Die achsenvariante Abbildung - The mapping with variant axis 

Die im Folgenden untersuchte Abbildung lasst sich 
nun nach folgender Vorschrift konstruieren: 

Gegeben ist der Kreis k(M,r) und ein Punkt P ^ M. 

1. Schritt: Bilde P durch Inversion auf den 

Hilfspunkt Q ab. 

2. Schritt: Bilde P durch Punktspiegelung an Q 

auf den Bildpunkt P’ ab. 

Man beachte, dass jeder Urpunkt durch diese 
Vorschrift sein „individuelles“ Spiegelungszentrum 
zugeordnet bekommt. 



Q is the image point of P applying the inversion wrt circle k(M,r). P* is the image point of P applying a 
reflection on Q. Each single point has its own reflection centre. 


Herleitung der Bildpunktkoordinaten - Deriving the coordinates of the image points 

Zur Herleitung einer Koordinatendarstellung fur den Bildpunkt P’ der Abbildung verwenden wir ein Koordi- 
natensystem mit dem Kreismittelpunkt als Ursprung. 

For finding the coordinates of P’ we use a system with the centre of the inversion circle as origin. The 
following calculation uses vector representation: 

Zunachst gelten die Vektorgleichungen: 


(1) MQ • MP = r 2 , da Q Bildpunkt von P bei der Inversion an k ist und (because of the inversion) 
— MP + MP' 

(2) MQ = , da Q als Spiegelungszentrum die Strecke PP’ halbiert. 



M 



Setzt man (2) in (1) ein, dann erhalt man: (*) MP • (MP + MP ’) = 2r 2 . Setzt man MP = 


und MP' = 



erhalt man zusatzlich zu dieser vektoriellen Abbildungsgleichung (**) 


x 

x' 


auf einer Geraden, die durch M, also den Ursprung geht. 


y_ 

y' 


, denn P und P’ liegen 


Einsetzen der Koordinaten in (*) liefert gemeinsam mit (**) ein Gleichungssystem fur die Bildkoordinaten x’ 
undyh DERIVE liefert als Bildpunkt von P(x,y): 


(Koordinatendarstellung der Bildpunkte) 


P' = 




U' 


2 r 2 x 
x 2 + y 2 

2 r 2 y 

2 2 

^x +T 


-x 


y 
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I- 1 1 gebralca lclotherlprgri I ole 1 Up! 1 

■[x y] ■+ mp : [xx yy] * rnp_ 

[xx yy] 

_ np + Dip- _ r 

x + xx y + yy 1 

■ 2 4 1 

2 2 J 

■ dotP^nq , np> = + equl 

2 

x ' 1 XX X 

y2 + = r 2 

2 2 

2 2 r 

■ x ■ yy = xx ■ y equ2 

yy ■ x = xx ■ y 


MMGH FiliD liUTD 


FiiF: 


Vj^Tnigebra 


Y F::t Y FHt ¥ F? T Ffi^ T "l 

i|Calc|Qiher|PrgnIQ|ciean Up| 1 


i x ■ yy = xx ■ y + equ2 yy ■ x = xx ■ y 

1 solue(equl and equ2,-Cxx yylO 
- 2 - 2 - y 2 ) x , (2 ■ r 2 - x 2 - y 2 ) ■ y 


and yy = 


.k: 


, 2 ) x 


x ^ + y^ 
= u and y = i 


MMGH FiliD flUTD 


FAR 19/30 


(You can find the DERIVE calculation in curve Jexicon.mth.) 


Konstruktion der Kurven - Construction of the Curves 


Zu den angekiindigten Kurven gelangt man nun, indem man Parallele zur y-Achse in verschiedenen Abstanden 
zum Ursprung punktweise abbildet. Bildet man eine Parallele zury-Achse mit Abstand u zum Ursprung ab, so 
haben die Geradenpunkte die Koordinaten (u,t) und werden demnach abgebildet auf Kurvenpunkte mit den 
Koordinaten: 


r 2 r 2 u 
K u 2 +t 2 


-u, 


2 r 2 t 
u 2 +t 2 


\ 

-t . 

J 


Mit t als Kurvenparameter liefert DERIVE als Bilder der Geraden - unter Verwendung von Schiebereglern fur 
den Radius r des Inversionskreises und den Abstand der senkrechten Geraden u : 




Die entstehenden Kurven ahneln nun sehr stark Kissoiden, Strophoiden und Trisektrices. Ob die abgebildeten 
Kurven aber wirklich die genannten sind lasst sich erst an Hand ihrer algebraischen Gleichungen entscheiden. 

We obtain the curves mentioned above as images of vertical lines x= u. The resulting curves are 
described by the parameter representation with parameter t. We can install slider bars for rand u and 
the resulting curves are reminding us on Cissoids, Strophoids and Trisectrices. Whether the 
conjecture is true or not can be decided by inspecting their respective algebraic equations ... 

But first of all let me - now in 2007 - reproduce the mapping using GeoGebra. The first figure shows 
the general construction. The we can edit the x-coordinate of point B which is the u from above and 
enter u = x = r ^2, u = x = r and u = x = r ! V2. 


_J Free objects 

■a A= (3.5, -2.03) 
■a B= (4.63,0) 

■a P = (4.63, 2.33) 

_3 Dependent objects 

-a M = (0, 0) 

-a P' = (i.oi, 0.51) 
a Q = (2.82, 1.42) 
a a: x = 4.63 
-a b = 4.06 
a i; x 2 +y 2 = 16.38 
O r= 4.05 
Auxiliary objects 



D-N-L#14 


Thomas Weth: A Lexicon of Curves (4) 


p!3 



x-u - r • V 2 



x = u = r 



Ca Free objects 

!■ Q 

i * 

[ 3 P = (4.05, 2.33) 


A =(3.5, -2.03) 

B = (r/sqrt(2), 


The Tl-Procedure: 


f » In 1 gebrata fclotherlprgn I die 1 ean Up! 1 

[2-n 2 - t 2 -u 2 )u 

+ xtl(t) Done 

" t 2 + u 2 

_ (2 ■ r 2 - x 2 - y 2 ) ■ y 

| x = u and y = i 

( 2 . r 2_ t 2_ u 2). t 

■M 

+ 

X 

1 


*2^,2 


MMQH RAD liUTD 

PAR 19/30 




ra n 


J-LQTS 

ytl: 


I z r - - t . - - u-j ■ t. 

^2 

^xt2=xtl(t) | r = 2 and u=|2-j2 

^yt2=yilCi) | r = 2 and u=|2-j2 
■^xi3= 2 ■ cost! tO 






f » 1 gebrata fclotherlprgn I die 1 ean Up! 1 

f 2 . r 2_ t 2_ u 2). t 

■ - = ? J + ytl(t) 

t, 2 + u 2 

■ soloe(x = xt-UtO and y = ytl(t)^ 

+ u" 

Done 

Cy 

.f u ”Ju (u + x) - 2 ■ r 2 j 

- 1 - -1 1 1 - ^ 

ana y — . 

■Ju + X 



MMQH RAD AUTD PAR 19/30 


f^fe^lT^ceW^aphlM^hlD^^IlV 1 1 
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Algebraische Kurvengleichungen - Algebraic Equations of the Curves 

Bezeichnet man der Einfachheit wegen die Koordinaten der Kurvenpunkte mit x und y so erhalt man die 


Parameterdarstellung 


2 r 2 u 


2 r 2 t 


x = - 


— t 


2 — = — 2 

^ U + 1 U + 1 j 


. Eliminiert man (mit DERIVE oder dem TI) den 


Kurvenparameter t , so erhalt man als algebraische Gleichung fur die Kurven: 


x 3 +x 2 


2 r 


2 A 


u -- 


+ y 2 (x + u2) — 0 (algebraic equation of the curves) 



Setzt man hierin speziell: 

♦ u = r-j2: so erhalt man als Kurvengleichung j/ 2 (r • V 2 + x) + x 3 = 0 , und nach einer Spiegelung der 

Kurve an der j-Achse (x — » -x): y 2 (rj 2 -x)-x 3 = 0, also die algebraische Gleichung der 

Kissoide. (Vgl. Lexicon of Curves (1)) 

♦ u = r: so erhalt man die Kurvengleichung: y 2 (r + x) + x 2 (x - r) - 0, also die algebraische Gleichung 

einer Strophoide. (Vgl. Lexicon of Curves (2)) 



so erhalt man als Kurvengleichung: y 2 



2 


-X 



= 0, was mit r = 1 die 


algebraische Gleichung einer Trisektrix von MacLaurin liefert. (Vgl. Lexicon of Curves 

( 3 )) 


Insgesamt liefert damit eine einzige Abbildung die aufgezahlten algebraischen Kurven als Bilder von speziellen 
Geraden. 


Das Problem: Bilder von Kreisen - The Problem: Images of Circles 

Die Untersuchung von Bildern von Kreisen unter der obigen Abbildung kann nun prinzipiell nach demselben 
Schema erfolgen. Allerdings ergeben sich hier algebraische Probleme, die ich trotz DERIVE nicht "in den 
Griff' bekam. Sicher bin ich, dass es sich bei den Bildern auch wieder um algebraische (und nicht etwa 
transzendente) Kurven dritter und vierter Ordnung handelt. Vermutlich fmden sich unter diesen Kurven die 
Lemniskate von Bernoulli, die bekannten Pascalschen Schnecken sowie als Spezialfall die Kardioide. Zum 
Abschluss seien noch einige Kurven dargestellt, die sich als Bilder von Kreisen mit der Parameterdarstellung 
(t + 1.5 cos a, 1.5 sin a) ergeben. Die Mittelpunkte der Kreise liegen also auf der x-Achse, die Radien sind 
r = 1.5, der Radius des Inversionskreises ist 1. 
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Interessante Kurven ergeben sich auch beim Abbilden von Kegelschnitten. Allerdings ist eine Herleitung der 
algebraischen Kurvengleichungen in diesem Fall sehr (zeit-) aufwandig. Ausfiihrlich ist die Untersuchung der 
oben angegebenen Abbildung und ihrer Eigenschaften dargestellt in "Didaktik der Mathematik, 2 (1991), 
S. 145 -164”. 

Allen, die sich mit dieser Konstruktion beschaftigen wollen, ware ich sehr dankbar, wenn Sie mich iiber ihre 
Entdeckungen und Ergebnisse informieren wiirden. 


Using this mapping we can investigate the images of circles. Here we are facing algebraic problems 
which I couldn't solve even using DERIVE. I am sure that the images are algebraic curves, too (of 
order three and four). I think that we will find the lemniscate of Bernoulli, the well known snails of 
Pascal with the Cardioide as a special case among them. In the figures you will find some curves 
which are images of circles with the parameter form [u + 1.5 cos(t), 1.5 sin(t)]. The centers of the 
circles are lying on the x-axis. The diameter of the inversion circle is 2. 

Interesting curves appear as are the results of mapping conic sections. Deriving the algebraic 
equations of the curves is a very time consuming work. You can find the mapping introduced above in 
"Didaktik der Mathematik, 2 (1991), p 145- 164". 

I would appreciate each investigation and result concerning this special mapping. Thank you! 


In the original version of 1994 Th. Weth added four figures, in 2007 I used the GeoGebra file from 
above and redefined the line x = u as the circle given above. In an instant I can see the new image 
pictures of the circles. 




The next page shows the DERIVE calculation followed by the graphic representation using a 
slider bar for varying parameter u. DERIVE and GeoGebra form a wonderful pair for 
investigating mappings and loci. 
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2 2 
x + y = 1 


3 3 

u + COS(t), SIN(t) 

2 2 


map := 


2-x 2-y 

x, y 


2 2 
Lx + y 


2 2 
x + y 


circ_img := 5UB5T(map, [x, y], circ) 


circ_img := 


2 2 2 
36 ■ u -COS(t) + 3.(12- u + l).COS(t) + 2-u.(4-u +1) 


2 ■ (12 ■ u -COS(t) + 4 ■ u + 9) 


3 ■ SIN(t) ■ (12 ■ u -COS(t) + 4 ■ u + 1) 


2 ■ (12 ■ u -C05(t) + 4 ■ u + 9) 


#14 : y = 


22 2 42 22 2 ~ 

(7(64 -u -x + 16-u-x-(4-u - 17) + 16- u -8-u +l)-8-u-x-4-u +1) -((8-u-x + 12 -u + l)-7(64- 


2 2 2 2 
128 -u ■ (7(64 ■ u -x + 16-u-x- (4-u - 17) - 

22 2 42 22 242 

u -x + 16-u-x- (4-u - 17) + 16-u - 8-u + 1) - 64-u -x + 128-u-x-(l - u ) - 80-u + 280-u - 1) 


4 2 2 2 

+ 16-u - 8-u + 1) - 8-u-x - 4-u + 17) 
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#20: 

ell := [u + a- COSCt), 

b-SIN(t)] 

#21: 

ell_img := SUBSTOnap, 

<v 
1 — 1 

i i 

#22: 

ell_img := - 



332 2 222 222 

a -03500 + 3 -a -u-COSOO + a-C0SOO-(b -SINCt) + 3-u - 2) + u-(b -SINCt) + u - 2) 


2 2 2 2 2 
a -COSOO + 2-a-u-03SCt) + b -SINCt) + u 

2 2 2 2 2 
b-SINCt) - (a -COSCt) + 2-a-u-QOSCt) + b -SINCt) + u - 2) 


2 2 2 2 2 
a -GOSCt) + 2-a-u-COSCt) + b -SINCt) + u 


The images of ellipses represented supported by slider bars in DERIVE and of general conics 
(GeoGebra): 
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Fluid Flow in DERIVE™ 

C.E. Reuther de Siqueira & J.V. Domingos, Petropolis, Brasil 


ABSTRACT 

The DERIVE program shows itself as a powerful tool to the teaching of Mathematics, 
Physics, Thermodynamics and Strength of Materials among others. The students of 
Engineering courses are hopeful with the results presented in DERIVE, mainly that with 
graphics visualization. In this contribution we use the graphic facilities of DERIVE for the 
graphic visualization of the fluid flow through a cylinder. 


INTRODUCTION 


The Navier-Stokes equations for viscous fluid flow is very complex for the mathematical analysis. 
The Computational Fluid Dynamics (CDF), is using Finite Difference Methods for the solution of 
these equations, together with methods of grid generation. But, simplification of Navier-Stokes 
equations could be done without a sensible loss of reality. So, one could find analytic solutions more 
easily. Here we will use one of these simplifications using the potential flow in polar coordinates and 
in steady incompressible flow. From this the continuity equation is [1,2] : 


V V = 


drV 8K 


d. + <30 


v 


= 0 


Figure (1) illustrates the cylinder studied and the speed of the fluid in infinity. 


( 1 ) 



Using the potential flow, the streamfunction, which satisfies continuity equation (1), could be [3] : 


y/ = U 0 r 


( 2 ^ 

1 a 

* 2 

V 

v ' y 


sinG 


( 2 ) 


where y/ ist he streamfunction and a 2 = with A being a constant and U () the initial speed of the 

^0 

fluid. The solution of equation (2), using DERIVE will be: 

¥ 


+ . 


r = 


U 0 sin0 ^ 


¥ 

U (} sin0 


+ 4a 


2 


( 3 ) 
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An idea of the streamlines could be obtained making y/= 2 a in equation (3) and the result will be: 

r = a |csc0 Wl + csc 2 0j (4) 

In equation (4) the function varies from infinite, in 0 = 0, to 1 in 0= tt 2, and then becomes ininity in 
0=0 and 0= n. The streamlines are presented in figure (2). 



Figure 2. The streamlines around the cylinder 


For figure (2) we use the following sequences in DERIVE: 


#1: 

Fluid Flow 


#2: 

Up lines 


#3: 

VECTORCa-CCSCCe) 

2 

+ /ci + cscce) 

#4: 

■For plotting: tt/1 

10 <- 6 <= 19n/2(M 

#5: 

Down lines 


#6: 

VECTOR (a -(C5C(e) 

2 

+ ./Cl + CSCC0) )), 

#7: 

The cylinder: 


#3: 

2.5 



CONCLUSION 


Equation (4) is written in polar coordinates. If DERIVE doesn't have the possibility for plotting these 
curves in these coordinates the problem would become even more complex. The development of 
numerical procedures and graphics for visualization of the streamlines would be necessary. In this 
contribution the speed and facility for visualization in the problem of Fluid Dynamics was shown. 
Like in this work others like Thermodynamics and Strength of Materials could be elaborated. Once 
more DERIVE showed itself as powerful tool. 
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I didn't translate or give a short English version of Mr Ramos' contribution, because I really believe 
that you will understand his ideas presented in his file. As a little support I changed his DERIVE- 
comments, giving short translations. Spanish is a wonderful language and I hope you will enjoy his 
desplazamientos. I really did. Muchas gracias! Josef 


Los desplazamientos 
en las funciones elementales 

P. Familiar Ramos, Valencia, Spain 

Este articulo es un estudio de las transformaciones geometricas en el piano aplicadas a una funcion 
u:=f{x), utilizando el programa DERIVE. Para ello se defmen las respectivas funciones que corresponden a los 
movimientos de traslacion, rotacion, reflexion, o reflexion en deslizamiento (las cuatro isometrias del piano ), y 
tambien se estudia la homotecia. 

A1 tratar de enfocar el problema, en seguida nos encontramos con una limitacion: dado que incluso 
una pequena rotacion aplicada a la funcion u=f(x) origina casi con seguridad que la funcion deja de ser 
uniforme... lo mejor es trabajar en parametricas al efectuar estas transformaciones. 

^Alguna otra limitacion? Si. Como DERIVE no es capaz de representar en parametricas funciones que 
tengan discontinuidades asintoticas, por ejemplo no dibuja (t,l/t) sabre un intervalo que contenga el cero (v.g.- 
3.1416 < t < 3.1416), nos limitaremos a considerar funciones continuas. 

Sin embargo, DERIVE es capaz de dibujar en parametricas funciones con discontinuidades de salto 
finito. For ejemplo representa sin problemas la funcion 

[x, IF(x < 0, 3x e x + 2, SIN(3x))] en -5 < x < 5, 



de modo que tambien podemos aplicar las transformaciones geometricas a este tipo de funciones. 

Es interesante experimentar con funciones sencillas y ver como actuan los movimientos. Presentamos 
ahora el listado que contiene las transformaciones geometricas, listas para ser usadas, y mostramos a 
continuacion algunos ejemplos concretos, que sin duda pueden ser mejorados, pero que sirven para ilustrar esta 
teoria. 

#+++##++++##+++##++++##+++###+++##+++###+++##++++##+++##++++##+++###+++##+++###+++##+ 

TRANSFORMACIONES GEOMETRICAS EN EL PLANO 

aplicadas a una funcion continua u:= f(x) 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

Traslacion de vector (a,b) - Translation by the vector (a.b) 

#1: TRASLACION Cu, a, b) := [x, u] + [a, b] 


Traslacion horizontal del longitud |a| - Horizontal translation 
#2: TRASLACION_H(u , a) := TRASLACION Cu, a, 0) 
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Traslacion vertical del longitud |b| - Vertical translation 
#3: TRASLACION_V(u , b) :i TRASLACION (u, 0, b) 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


Rotacion de centra (xO.yO) y amplitud a (en grados) - 
Rotation with rotation center (xO, yO) by a degrees 


#4: 


RQTACIQN(u, xO, yO s «) := [x - xO, u - yO]- 


COS&x-l 0 ) SIN(qe- 1°) 
- SIN(ot -1°) CQS(of-l 0 ) 


+ [xO, yO] 


Simetria central respecto del punto (xO.yO) - Reflection with center (xO.yO) 
#5: SEMIGIRO_(u , xO, yO) := RGTACIONCu, xO, yO, ISO) 

Giro de amplitud a alrededor del origen - Rotation by a deg around origin 
#6: GIROCu, of) := RGTACIGN(u, 0, 0, of) 

Simetria central respecto del origen - Reflection with origin as center 
#7: SEMIGIRO(u) := GIRG(u, ISO) 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


Reflexion respecto de la recta ax + by + c = 0 - Reflection with respect to the line ax + by + c = 0 

2 - 3 .- b 


#S: REFLEXIONCu, a, b, c) := [x, u]- 


2 2 
a - b 


2 2 
a + b 


2 - a - b 

2 2 
a + b 


2 2 
a + b 

2 2 

a - b 

2 2 

a + b 


Simetria respecto del eje Ox - Reflection with respect to the x-axis 
#9: S IM ETR I A_OX ( u ) := RE FLEXION (u, 0, 1, 0) 

Simetria respecto del eje Oy- Reflection with respect to the y-axis 
#10: SIMETRIA_GY(u) := RE FLEXION (u , 1, 0, 0) 

Simetria respecto de la diagonal principal - Reflection with respect to y = x 
#11: S IM ETR I A_DI AGO N AL ( u ) := RE FLEXION Cu, -1, 1, 0) 

Simetria respecto de la diagonal secundaria - Reflection with respect to y = -x 
#12: S IM ETR I A_DI AGO N AL 2 ( u ) := RE FLEXION (u, 1, 1, 0) 


2-a-c 

2 2 
a + b 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


2-b-c 


2 2 
a + b 


Homotecia de razon k con centra en el punto (xO.yO) - Enlargement by factor k with center (xO.yO) 
#13: HQMOTECIA_Cu, xO, yO, k) := k - [x - xO, u - yO] + [xO, yO] 

Homotecia de razon k centrada en el origen - Enlargement by factor k with center (0,0) 

#14: HQMOTECIACu, k) := HOMOTECIA_Cu , 0, 0, k) 


*+*+*+*++++++++++*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+ 
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Reflexion en deslizamiento del longitud |d| respecto de la recta ax+by+c=0- 
Reflection and glide by |d| with respect to the line ax + by + c = 0 

#15: de la recta ax + bx + c = 0 

b-d 

DESLIZAMIENTO (u, a, b, c, d) := RE FL EXION (u, a, b 5 c) + 

#16: 22 

_ 7(a + b 3 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


a-d 


2 2 

■/(a + b ) . 


Ejemplos - Examples 

The original curves are black and the results of the transformations are red. 


2 

#17: - x 

u := e 


#13: 


#19: 


TRASLACIGNCu, 2, -1.53 

T 2 1 


- x 

x + 2, e 


3 

2 


# 20 : 


r 2 

- x 


TABLETS , x, -5, 5, 0.01J 
#21: (TABLE (TRASLACION (u , 2, -1.53, x, -5, 5, 0.0133-U2 



u := 2-SIN(x3 - COS(2-x3 


S IM ETR I A_OX ( u 3 


GIRG(u, -303 
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Reflection 


Reflection 



SIMETRIA_DIAGONALCu) 


Reflection & Glide 



u := 


3- X 


1 + X 


0 ■ x + 1-y - = 0 

2 


DESLIZAMIENTG 


u, 0, 1, - 


1 > 

, 2 

2 j 



3- x 


u :i 


1 + x 


2 ■ x - 3-y - 1 = 0 
DESLIZAMIENTOCu, 2, -3, -1, -2) 



2 

u := SINCx ) 



2 

u := SINCx ) 



^ i ^ 


r 3 '' 

□= 
o 
zc 
o 
— 1 
m 
n 

H 

P 3 - 

u, 0, -4, 

□= 
o 
zc 
o 
— 1 
m 

n 

H 

P 3 - 

u, 2, 2, 


V 


J 
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Ahora que conocemos major las notaciones, podemos ir mas alia jugar un poco con expresiones del 
tipo 


VECTOR (GIRO (SIN (x) , k) , k, 11.25, 180, 11.25) 
obteniendo en este caso: 



^Qien iba a decir que este resultado es un «conjunto de senos»l 

It is not so difficult to transfer the transformations on the CAS-TI devices. However there are some 
specialities which must be considered. I convert the matrix into a list - then it is easier to store the two 
components of the parameter representation in the Y=EDITOR. Piecewise defined functions make 
problems by transfering them into the Y=EDITOR (see below). 


matDl i st ( [t , u] + [a , b] ) | x=t -> trasl (u , a , b) 


UpI ] 

PiSfeafeSitrWi W5i.t s, .1 1 


■ natP 1 ist([t u] + [a b]) | x = t + trasl (u , > 

Done 

■ trasl U " x2 , 2 , -l.s) 

{t + 2 e" t2 -3x2 } 

■{t + 2 e " t2 - 3^}[l] + xt2(t) : {t + 2 ► 

Done 

J'LDTl" 
r rue 1 ■ 

^xtl=t 

-t 2 

^ytl=e 2 

■^xt2={t + 2 - 3x2/ . 1. 

v 'yt2={ t + 2 e _t - 3 x 2 ) [ 2 ] 
xt3^ 







...Ct+2.e A <-t A 2>-3/2>[2]^t2<t> 

xt3<t >= 

MAIN RAD EKACT PAR 3/30 

MAIN RAD EKACT PAR 

MAIN RAD EKACT PAR 


matDlist( [x-x0 , u-y0] * [cos (a) ,sin(a) ; -sin (a) ,cos(a) ]+[x0,y0] ) |x=t -» rotac(u,x0,y0,a) 


|^lu«:VlloCC^iob;, ' .1 ] 

|t 1 gebrafca fclotherfprgn I ole 1 ean Upf 1 


Done 

■ rotac(2 ■ sin(x) - cost 2 ■ x), 1, 1,60°) 
pcp| ( 2t)_^. £in(t)+ t + ^ + ^ 2 

.poo,(20_^. s . naJ+ t + ^ +w * 

* L 'sin(ot) costal 

Done 

■ rotac(2 ■ sintx) - cost 2 ■ x), 1, 1,60°) 
pcos(2«.^. £ . nW+ tp +1 ^ .> 

Done 



vG] ] > 1 x=t->rotac<u,xG,vG,a) 

...Ml]+xt4<0:ans<lM2]-Wt4<0 

MAIN RAD EKACT PAR 3/ 6 

MAIN RAD EKACT PAR fi/30 

MAIN RAD EKACT PAR 
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matDl i st ( [x , u] * [ - (a A 2-b A 2) , - 2a*b ; -2a*b , a A 2-b A 2] / (a A 2+b A 2) + 
[ -2a*c , 2b*c] / (a A 2+b A 2) ) | x=t reflex(u,a,b,c) 



Substituting t for x does not apply for the piecewise defined function. I helped myself by defining the 
STEP-function, which returns 0 for x<0 and 1 for x>0. You could also define u as a function of t. 
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See finally the family of rotated sine waves on the Tl-screen. Plotting takes some time! 


■^xtl = ^ cos( k) ■ t - sin( k) ■ sin(t) cost; k) ■ sii> 



vtl <t >=<cqs <k>^t~sin <k>^sin<t„, 

MHIH FiliD flUTD PmF; 



In my opinion it might be a nice task for students to transfer the transformations on the Tl trying to 
overcome all problems. There are transformations left (eg shear). So they could also accomplish the 
collection of planar transformations. Josef 


Aufgaben aus der Elektrotechnik - nicht nur 
fur den Unterricht in beruflichen Schulformen 
Problems from electrical engineering - not only for Vocational Schools 

H. Scheuermann & N. Wild, Hofheim, Germany 


Mit diesem Beitrag mochten wir anwendungsorientierte Beispielaufgaben aus der Elektrotechnik 
fur den Schulunterricht unter besonderer Beriicksichtigung von DERIVE™ vorstellen. Die 
Aufgaben stammen aus unseren Unterrichtskonzepten fur die beruflichen Schulformen 
Berufliches Gymnasium und Fachoberschule mit den Schwerpunkten Elektrotechnik bzw. 
Informatik sowie (nach angemessener didaktischer Reduktion) fur die Teilzeitschulform 
Berufsschule (insbesondere Kommunikationselektronik). Alle Beispiele wurden von uns im 
Unterricht dieser Schulformen meist mehrfach ausprobiert und evaluiert. Aufgrund der - nach 
unserer - Einschatzung positiven Unterrichtserfahrungen glauben wir, dass sich die folgenden 
Anwendungsaufgaben - zumindest partiell - auch als Unterrichtssequenzen fur den Mathematik- 
oder Physikunterricht an den allgemeinbildenden Schulen eignen. 1 

1. Aufnahme der Kennlinien eines belasteten Spannungsteilers 2 
Characteristics of a loaded voltage divider 


A ufgabenstellung: 

Die Kennlinien eines belasteten Spannungsteilers 
(vgl Abb. 1) sollen graphisch dargestellt werden. 
Durch Veranderung des Lastwiderstands R L soil 
dessen Einfluss auf den Verlauf der Kennlinien 
untersuchtwerden. 

Problem: Give the graphic presentation of the 
characteristics of a loaded voltage divider. 
Investigate the influence of the load resistor R L on 
the characteristics. 


Ri 


U 


R 2 


Hl 


Rl 


Abb. 1: Belasteter Spannungsteiler 


Das DERIVE-Listing 3 konnte folgendermaBen aussehen 4 - bis zur normierten Form (Ausdruck #12) 
The DERIVE-procedure could look like as follows - up to the normalized form (expr # 12) 
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Calculation of the total resistance Rg 

with Rp = R1 + R2 as resistance of the potentiometer. 

#1: CaseMode := Sensitive 

#2: InputMode := Word 


R2-RL 

#3: Rg := R1 + 

R2 + RL 

R2-RL 

#4: Rg := CRp - R2) + 

R2 + RL 


U 

#5: Ig := 

Rg 

U-(R2 + RL) 

Ig := 

#6: 2 

R2 - R2 - Rp - RL-Rp 


Ig-R2-RL 

#7: UL := 

R2 + RL 


R2-RL-U 

UL := 

#3: 2 

R2 - R2 ■ Rp - RL-Rp 

Division by U and normation by substitution: 
F(x)=UL/U, R2 = x* Rp and RL = Rp/z 

UL R2-RL 


#9: U 2 

R2 - R2 - Rp - RL-Rp 

Rp 

tx-Rp) 

UL z 

# 10 : = 

U 2 Rp 

(x-Rp) - (x - Rp) - Rp Rp 


UL x 


#11 : U 2 

X -z - x-z - 1 

Domain: [0,1] 

We plot some characteristics by variation ofz 

and now with DERIVE 6 we add a slider bar to animate the characteristic! 


F(x , z) := 


#12: 

If 0 < x < 1 



- x/(x A 2-z - x-z - 1) 


#13: 

VECTOR(F(x, z), z, [0, 0.8, 

2, 4, 8, 16, 40]) 


TABLE ''FCx, z) := 

, x, 0, 1, 0.001 A 

#14: 

If 0 < x < 1 



- x/(x A 2-z - x-z - 1) 



Abb. 2: Kennlinien des belasteten Spannungsteilers 

Durch Variation von z (hier z e {0, 0.8, 2, 4, 8, 16, 40}) erhalt man die Kennlinienschar im Plot- 
Fenster (vgl. Abb. 2). 

Variation of z - z e {0, 0.8, 2, 4, 8, 1 6, 40} - gives the family of characteristics. 
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2. Die belastete Spannungsquelle - The loaded supply point 


A ufgabenstell ung: 

Gegeben ist eine Spannungsquelle mit dem In- 
nenwiderstand Rj und einem regelbaren Lastwi- 
derstand R L (vgl. Abb. 3). Die Leistung P L soil 
in Abhangigkeit vom Lastwiderstand Rl 
untersucht werden. 

Problem: Given is a supply point with an internal 
resistance R| and a regulable load resistor R L . 
Investigate the power P L in dependence on R L . (It 
might be easier for the students to work with 
special values for the source voltage U 0 and for 
R.)- 


i 



Wie bei den folgenden Aufgaben ist auch bei dieser die Wahl eines geeigneten, adressatengerechten 
Anspmchsniveaus unabdingbar. In diesem Sinne konnten die inhaltlichen Anforderungen reduziert 
werden, wenn fur die Quellenspannung und den Innenwiderstand bereits in der Aufgabenstellung 
Zahlenwerte angegeben werden. Diese didaktische Reduktion ermoglicht die Koordinaten fur die 
maximale Leistung ohne groGeren rechnerischen Aufwand aus dem Funktionsgraphen mit Hilfe des 
Cursors auszumessen. Dieser Vorgehensweise entspricht u.a. der Losungsvorschlag in Abb.4 5 . 

Die Auswertung der Zeilen #1 bis #9 (inkl. Funktionsgraph) legen die Annahme nahe, dass die Lei- 
stung am Lastwiderstand maximal wird, wenn R L = Ri gilt. Die Verallgemeinerung und damit die 
Richtigkeit der Annahme beweisen die folgenden Zeilen im Listing (vgl. Abb. 4 5 ). 

Interpretation of expressions #1 to #9 including the graph gives the idea that power will be maximal at 
the load resistor. Generalization and with it the truth of this assumption can be shown in the following 
listing -see fig. 4 5 . 

Eine alternative Losungsvariante weist das folgende DERIVE-Listing auf. Da hier auf das Instrumen- 
tarium der Differentialrechnung verzichtet werden kann, eignet sich dieser Losungsweg auch fur die 
llte Klasse des BGs in der Oberstufe. Diesem Losungsansatz liegen pragmatische Uberlegungen 
zugrunde. Die Leistung P L am Lastwiderstand kann als Produkt aus Spannung und Strom beschrieben 
werden, also P L = U L I. Die Variation des Lastwiderstands Ri e [0, go] bewirkt eine Anderung beider 
Faktoren. Die Grenzwerte von R L zeichnen sich durch ihre praktische Bedeutung aus: 

(I) R l = 0Q (Kurzschluss) — » der Strom I wird maximal, die Spannung U L = 0V — » P L = 0W 

short circuit 

(II) RL — » go („offene“ Klemmen) — » die Spannung wird maximal, der Strom I = 0A —> PL = 0W 

open binders 

Da fur 0 < P L < oo die Leistung P L > 0 ist, muss mindestens ein R L existieren, das fur P L ein Maxi- 
mum besitzt (dies ist anschaulich klar; der mathematische Hintergrund, dh der Satz von Rolle, sollte 
aber in diesem Zusammenhang kein Thema fur den Unterricht sein). Das Auflosen der Gleichung 
nach R l in Zeile #3 fiihrt zur Losung einer quadratischen Gleichung mit positiven bzw negativen 
Wurzeltermen (Zeile #4). Die Existenz von zwei Losungen wird in der graphischen Darstellungs- 
bene plausibel, wenn das Diagramm in Abb. 4 durch eine Parallele zur ersten Achse 
(mit 0 <P L <P l max) erganzt wird. (Die Sekante weist mit R L i und R L 2 zwei Schnittpunkte auf). Durch 
Verschieben dieser Parallelen in positiver Richtung der zweiten Achse geht die Sekante in eine 
Tangente (mit P L = PLmax und dem Argument R L i = Rl 2 ) iiber. Bezogen auf die Losungen der 
quadratischen Gleichung bedeutet dies, dass die Diskriminante gleich Null wird (Zeilen #5 ff). 
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#1: [CaseMode := Sensitive, InputMode := Word] 


#2: PL = 


UO 


RI + RL 


-RL 


#3 : SOLVE 


PL = 


UO 


RI + RL 


-RL, RL 


#4: RL = - 


2 2 
JW - 4-PL-RI) • | UO | 2-PL-RI - UO 


v RL z 


7(U0 - 4-PL-RI)- | UO | 2-PL-RI - UO 


2 -PL 


2 -PL 


2 -PL 


2 -PL 


Set the discriminant = 0 


# 5 : 


2 

SOLVE (UO - 4-PLmax-RI = 0, PLmax) 


2 

UO 

PLmax = 

4-RI J 


expr #3 = expr #6 and solve for RL: 


#6: SOLVE 


UO 


U RI + RL 


UO 


•RL 


4-RI 


RL 


v UO = 0) 


Compare the DERIVE screens from 2007 and 1994: 
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3. Der Effektivwert eines sinusformigen Wechselstroms - 
The effective value of a sinusoidal alternating current 

A ufgabenstell ung: 

Beweisen Sie, dass der Effektivwert 6 eines sinusformigen Wechselstroms berechnet werden kann 
durch die Formel 

i ^max 

e// ' V2 • 

Problem: Proof that the effective value of a sinusoidal A.C. can be calculated by the formula given 
above. 


#1: [InputMode := Word, CaseMode := Sensitive] 

Labour of direct current at an ohmic resistor R during a period T: 

#2: [w_ := P_-T, P_ := R-Ieff ] 


2 

#3: W_ = Ieff -R-T 

Function for alternating current 

#4: I(t) := Imax-SIN(w-t + 4 ) 

Labour of A.C. at an ohmic resistor R during a period T: 


# 5 : 


Ww 


T 2 

I Pw dt, Pcjo := R-I(t) 

J 0 


The definition of the effective value leads to equation 
#6: W_ = Wbo 

2 2 2 
2 Imax -R-SINC2-T-W + 2-cf>) Imax -R-SIN(2-4) Imax -R-T 

#7: Ieff -R-T = + + 

4-w 4-w 2 


Substitute the angular frequency by 2 jt/t 


Imax -R-SIN 


2-71 

2-T + 2-ct? 

T 


#3: Ieff -R-T = - 


2 2 
Imax -R-SINC2-4) Imax -R-T 
+ + 


2-jt 


2-jt 


2 Imax -R-T 

#3: Ieff -R-T = 


#10: SOLVE 


2 ^ 
2 Imax -R-T 

Ieff -R-T = , Ieff 


■J 2 - Imax 


#11: Ieff = - 


2 


v 


2 


yR = 0vT=0 




p32 


H. Scheuermann & N. Wild: Elektrotechnik 


( 1 ) 


D-N-L#14 



Abb. 5: Wechselstrom AC = IQ(t), Effektivstrom I_(t) 


4. Einstellung des optimalen Arbeitspunkts bei symmetrisch angesteuerten 
Transistorstufen - Adjustment of the optimal operating point at symmetric steered 
transistor steps 


A ufgabenstell ung: 

Gegeben ist eine Transistorstufe laut Schaltbild (vgl. 

Abb. 6). Vom Transistor sind das Ausgangskennlinien- 

feld (I c = f(U C E)j zB als Derive-Listing) und die maxi- 
male Verlustleistung P = 45 W bekannt. 

(I) Bestimmen Sie den Arbeitspunkt der Transistor- 
stufe fur eine maximale Leistungsverstarkung mit 
U B , Rc als Parameter (Ansteuerung mit einem 
symm. Wechselspannungssignal). 

(II) Die Betriebsspannung U B betrage 12V. Berechnen 
Sie R c unter Beriicksichtigung von (I). 

(III) Zeichnen Sie den Graph fur die Widerstandskenn- 
linie in das DERIVE-Ausgangs-kennlinienfeld und 
iiberprufen Sie die Ergebnisse aus (I). 



Abb. 6: Transistorschaltung 


Problem: Given is a transistor step (fig. 6). We know the family of output characteristics of the 
transistor and its maximum power loss P = 45W. 

(I) Find the operating point for a maximum power amplification. 

(II) Operating voltage = 12V. Calculate R c under consideration of (I). 

(III) Plot the characteristic of the resistor to the family of output-characteristics and check the results 
of (I). 


#1: [CaseMode := Sensitive, InputMode := Word] 

Task (I) 

Equation of the hyperbola for the maximum powerloss 
P 


#2: Ip(Uce) := 


Uce 
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The 1st derivative of Ip(Uce) is Ipl 


#3: Ipl(Uce) := Ip(Uce) 

d Uce 


Ipl(Uce) := 


Equation of the resistor line 

1 Ub 

#5: Iw(Uce) := Uce + 

Rc Rc 

1st derivative of Iw is Iwl 


#6: Iwl (Uce) := Iw(Uce) 

d Uce 


#7: Iwl (Uce) := - 


The coordinate Uce in the osculation poii 


#3: [Ip(Uce) = Iw(Uce) , Ipl(Uce) = Iwl(Uce)] 


Ub - Uce P 


#9: Uce 


2 Rc 


Solve this system for Uce (= volatge in the operating point) and Rc 


Ub - Uce P 


L [Uce, Rc] 


2 Rc 


Ub Ub 2 

#11: Uce = a Rc = a Rc-Uce t 0 a Rc-Uce t 0 

2 4 - P 

Substitute for Uce in #2 = (collector) current Ic in the operation point 


Ub T 2 - P 


#12: Ip 


Task (II) 

Replace Uce in #9 (1st component) by Ub/2, solve for Rc and substitute the given values 
P Ub - Uce 


2 - P Ub 


Ub 2 - Rc 
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#16: 


#17: 


SOLVE 


2 - P 


Ub 


2 

Ub 

Rc = 

4 - P 


2 

12 

#13: r c = 

4-45 


#19: Rc = 0.3 


Ub 

, Rc 

2-Rc 


Task (III) 

Resistorline: replace Rc in #5 by Rc of #1 7 
Ub - Uce 

#20: Iw(Uce) = 

Rc 

Ub - Uce 

Iw(Uce) = 

2 

#21: Ub 


4 - P 

Ub - Uce 4-P-CUb - Uce) 


#22: Rc 2 

Ub 

4-45 -(12 - Uce) 


#23 : 2 

12 



Abb. 7: DERIVE™-Graphik: Ausgangskennlinienfeld. Verlusthyperbel. Widerstandsgerade. I c = 1(1)^) 


Original screen shot from DNL#14 


1 Den fachlichen Hintergrund beschreiben zB [1] und [2]. 

2 Eine ausfiihrliche kommentierte Unterrichtssequenz iiber den belasteten Spannungsteiler in [6]. 

3 Die Bezeichnung DERIVE-Listing beinhaltet die Gesamtheit aller Schreibzeilen zur Losung einer Aufgaben- 
stellung. 

4 Selbstverstandlich konnen bzw. miissen bereits bei der Aufgabenstellunge der Lerngruppe entsprechende 
angemessene Vereinfachungen wie zB konkrete Zahlenwerte fur U, R P , R L vorgenommen werden. 

5 Das Cursor-Kreuz im Plot-Fenster weist das Maximum des Graphen der Funktionsgleichung P L =f(R L ) aus; die 
Koordinaten des Kreuzes konnen links unten abgelesen werden, also x = R LM » 0,8 und y = Pl(Rlm) ~ 45W. 

6 Effektivwert: Der Effektivwert I eff eines Wechselstroms ist der sich aus den Augenblickswerten I(t) ergebende 
Dauerwert, der in einem ohmschen Widerstand die gleiche Warmearbeit erzeugt, wie ein Gleichstrom der 
gleichen Hohe (aus [4], Kapitel 2, S 33). 


will be continued in DNL#15 
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Josef Lechner realized J.H. Conway's famous Game of Life for Derive. I copy the first lines of the 1994 
DERIVE code together with one result. In the recent Derive version we can use the Derive Unicode 
font, which shows the “living cells” in a better way. Josef 

Conway 's CAME of LIFE 

Josef Lechner, Viehdorf, Austria 

(In 1994 the SUB-Operator was not available. We had to enter ELEMENT(a, i , j) instead of a 4,1 rj .) 


#1: LIFE - A Game from J.H. Conway 

#2: An attempt with DERIVE (Josef Lechner) 

#3: We have three creatures in our miniworld 



#5: 


CREA(a, i, j) : 
If aii-l-J = 1 
"Yes 11 
M No M 


#6: CREA(wor Id , 2,3)= Yes 


#7: CREA(wor Id , 1, 1) = No 


' 0 0 0 0 0 

0 0 10 0 
#31: t := 0 1110 

0 0 10 0 

. 0 0 0 0 0 

#32: EV0L(t , S) 



' 0 

0 

0 

0 

0 ' 


' 0 0 

0 

0 

0 ' 


' 0 

0 

1 

0 

0 ' 


' 0 

0 

1 

0 

0 ' 


' 0 

1 

1 

1 

0 ' 


' 0 

1 

1 

1 

0 ' 



0 

0 

1 

0 

0 


0 1 

1 

1 

0 


0 

1 

0 

1 

0 


0 

1 

1 

1 

0 


1 

0 

0 

0 

1 


1 

0 

1 

0 

1 



0 

1 

1 

1 

0 

, 

0 1 

0 

1 

0 

, 

1 

0 

0 

0 

1 

, 

1 

1 

0 

1 

1 

I 

1 

0 

0 

0 

1 

, 

1 

1 

0 

1 

1 

, 


0 

0 

1 

0 

0 


0 1 

1 

1 

0 


0 

1 

0 

1 

0 


0 

1 

1 

1 

0 


1 

0 

0 

0 

1 


1 

0 

1 

0 

1 



. 0 

0 

0 

0 

0 . 


. 0 0 

0 

0 

0 . 


. 0 

0 

1 

0 

0 . 


_ 0 

0 

1 

0 

0 . 


. 0 

1 

1 

1 

0 . 


. 0 

1 

1 

1 

0 . 



0 

1 

1 

1 

0 1 

' 0 

1 

1 

1 

0 ' 


' 0 

1 

1 

1 

0 " 

1 

0 

0 

0 

1 

1 

0 

1 

0 

1 


1 

0 

0 

0 

1 

1 

0 

0 

0 

1 , 

1 

1 

0 

1 

1 

i 

1 

0 

0 

0 

1 

1 

0 

0 

0 

1 

1 

0 

1 

0 

1 


1 

0 

0 

0 

1 

0 

1 

1 

1 

0 J 

. 0 

1 

1 

1 

0 . 


. 0 

1 

1 

1 

0 . 

took 

114. 8 

i se< 

i in 1994, 

ta 

kes 

now 0. 

53 

sec 

1 





You can find numerous great websites dealing with the Game of Life: 


http : // www . radicaleye . com/lif epage/ 
http : //www . ibiblio . org/lif epat terns/ 

http://de.wikipedia.org/wiki/Conways Spiel des Lebens 
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There are three functions for the game: 

LIFEfa, i, j) := 

If num(a, i, j) = 1 a (N(a, i, j) = 2 v N(a, i, j) = 3) 
# 11 : 1 

f 

DEATH (a, i, j) := 

If num(a, i, j) = 1 a (N(a, i, j) > 3 v N(a, i, j) < 2) 
# 12 : 1 

f 

BIRTHfa, i, j) := 

If num(a, i, j) = 0 a N(a, i, j) = 3 
#13: 1 

f 


NG(a,i j) determins whether the cell is living in the next generation or not. 

NG(a, i, j) :: 

If num(a, i, j) = 1 a (N(a, i, j) = 2 v N(a, i, j) = 3) 
1 

#14: If num(a, i, j) = 0 a N(a, i, j) = 3 

1 
f 


NG1(a) finds the full next generation. 

#15: NGl(a) := VECTOR (VECTOR (NG( a, 


#16: 


Example 1: 


#17: 


f 1 f 
f 1 f 
L f 1 f 


-i), -i, 1, 3), j, 1, 3) 


#1S: 


#19: NGl(s) = 


We want to observe the evolution of the population of miniworld s 

NG2(a) := VECTOR (VECTOR (NG (a, j, i), i, 1, DIMENSION(a )) , j, 1, DIMENSION(a)) 
# 20 : 1 

#21: EVOL(world, steps) := ITERATE5(NG2(m) , m, world, steps) 


#22: EVOL(s , 3) = 
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Example 2: 


#23: t := 


f f f f f ' 

f f 1 f f 

f 1 1 1 f 

f f f f f 

f f f f f . 


#24: EVDL(t, 8) 




This took 115 seconds in 1994, now it needs only one second for evaluation. A famous pattern is the 
following - see first one gliding step and then starting the GLIDER from left above gliding gently down: 


This is the famous GLIDER 

f f f f f f 

f f i f f f 

f f f i f f 

f i i i f f 

f f f f f f 

f f f f f f 


#26: gl := 


#27: EVOLtgl, 4) = 


I allocate the glider in the left upper corner of our 6 by 6 miniworld and then start the gliding process: 
You can experiment with various patterns starting in a larger world - see the “globe”. 


#29: g!2 


f 1 f f f f 
f f 1 f f f 
1 1 1 f f f 
f f f f f f 
f f f f f f 
f f f f f f 


#30: 


|EV0L(gl2, 12)1 
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#52: globe 


fffffffffff 
fffffffffff 
fffffffffff 
fffffffffff 
fffflffffff 
ffffl 1 lffff 
fffflffffff 
fffffffffff 
fffffffffff 
fffffffffff 
fffffffffff 


#55: IdMMSIBigJMAI 



As Josef Lechner wrote in a private communication he used Conway's Game of Life in Information 
Technology as an example for programming with two-dimensioned arrays. 

I remembered that I should have a paper on the Game of Life among my many collected articles from 
various newspapers and journals and finally I found an article from Scientific American 1984 where 
another cellular automat was presented: 


Josef Lechner: CONWAY's Game of Life 
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The Fredkin-Automat: 


Only the neighbour cells in horizontal and vertical direction are responsible for the next generation. 
If the number of "living" neighbour cells is even then the actual cell will die, 
if the number is odd, then the cell remains living or comes into life. 


#35: 

NF(a, i , j) := NB(a, 

i - 1, j) + NB(a, 

#36: 

NFG(a, i , j) := 

If H0D(NF(a, i , j; 
L 
F 

), 2) = 1 

#37: 

NFG2(a, i, j) := VECTOR (VECTOR (NFG (a, j 

#38 : 

EVDLF (wo rid , steps) 

:= ITERATES(NFG2(ra) 

#39: 

1EV0LF (globe, 



#41: EVOLFCgl, 4) = 


Have a 21 by 21 globe with two living cells 
and follow how this world will develop. 


The two Josefs wish much fun creating your 
own patterns of life. 


#57: gtobe_t 2 := 


http://www.univ-lemans.fr/enseignements/physique/02/divers/fredkin.html 
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Titbits from Algebra and Number Theory (2) 

by Johann Wiesenbauer, Vienna, Austria 

Hi, out there! Josef Bohm did a fabulous job when acting as my ghost writer last time, but here I am 
myself -eius ipsissimus, as it were. Before going to new stuff, I'd like to add a few comments to the 
DERIVE listing in DNL #13 (revised version) on page 31. Its main purpose was to show that Euclid's 
algorithm for the computation of the greatest common divisor of two natural numbers a and b (in 
symbols gcd (a,b)) can be easily implemented in DERIVE - even in its extended form which is often 
used to solve linear congruences, as shown below. 

For didactic reasons I started with the function euclid(a,b) which yields the divisions of Euclid's 
algorithm for a and b in a neatly arranged form. In (5) you can see its output for the values 
a = = 1234567891 and b = 234567891. In particular, the last nonvanishing remainder, which is 1 in 
this case, is the gcd (a,b). (Here and hereafter, if necessary , consult any textbook on Number Theory 
for the mathematical background, for example [1].) 


What about the number of steps needed for 


eucl id(fibl(14) . fibl(13 


Euclid's algorithm? It is easy to see that it is 

' 377 = 1 

233 + 144 ' 

strongly connected with the magnitude of the 

233 = 1 

144 + 89 

quotients which are all 1 in the worst case as e.g. 

144 = 1 

89 + 55 

in (19). Therefore you get an upper bound for the 

89 =1 

55 +34 

number of steps needed to calculate the gcd (a,b) 

55 =1 

1.A — 1 

34 + 21 

71 ± 

for a,b < s by counting the number of steps of the 

34 — -L 

21 =1 

+ -L7 

13 + 8 

algorithm of Euklid for the largest adjacent pair 

13 = 1 

8 + 5 

of the Fibonacci sequence 

8=1 

5 + 3 

0,1,1,2,3,5,8,13,21,34,... 

5 = 1 

3 + 2 

below s. 

3=1 

. 2=2 

2 + 1 
1 


Exactly this is done in (9) - (26) for 5 = lO 100 . (Note that the disappointing calculation time in (26) is 
by no means representative for actual applications of Euclid's algorithm since it is caused almost 
totally by the unnecessary storing of the whole tableau! In 1994 it needed 65 sec, now we need only 
0.05 sec.) In the following the result is derived in another way. 

For this purpose the explicit formula for Fibonacci numbers is needed, namely 

„ 1 rri+VJY fi-VsYl 


m 2 


1^1. 0.618 <1 andfi^l ,^5 
2 2 


this yields 


F„ * 


2 
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On the other hand, n - 2 is exactly the number of steps for the calculation of the gcd(F„ ,F n _ i ) (see 

above). Therefore its value is about Jog s , where X = is the number known as the ’’golden 

ratio” (cp. (24), (26)). In particular, it depends only linearly on the number of digits of s which is the 
best we could expect! 

On the other hand, the number of steps you get on average are not much better than those in the worst 
case (see below). It can be shown that in some sense the most likely value is about 


12 In 2 t 

2 — * ns 

71 


(again this n\) which is about 40% of the maximum. 


100 100 


DIM (eucl -id (RANDOM (10 ) , RANDOM (10 ))) = ISO 

100 100 


VECTOR (DIM (euc lid (RANDOM (10 ) , RANDOM (10 

12 
2 

k. 7T 


12 ■ LN(2) 100 

LN(10 ) 


100 


1 LN(2) 


))), k, 5) = [198, 169, 17S, 205, 191] 
- 100 ■ LN(5) 


12 ■ LN(2) 100 

LN(10 ) = iKEUBEWa^ 

2 

7T 

As for the EEA applied to a and b , each line 

fo-i n X; y.^ y { ] 

of the output (cp. 31)) has the following meaning 

r x - x { _ x a + y { _ x b , r { = x x a + y { b. 


In the last row, apart from r n = gcd (a,b) (2. entry) also the coefficients x n and y n (4. and 6. entry, 
respectively) can be seen. Again, the whole tableau is presented here only for didactical reasons, 
usually we are only interested in a subset of these three numbers! For example, we could use the EEA 
to compute the inverse of a mod m, in case gcd(a,m)=l: 

INV(a,m) == IF(GCD(a,m) = 1, 

MOD(ELEMENT(ITERATE(IF(MOD(ELEMENT(x,l), ELEMENT(x,2)) = 0, x, 
[ELEMENT (x,2), MOD(ELEMENT(x,l), ELEMENT (x,2)), ELEMENT (x,4), 
ELEMENT(x,3) - FLOOR(ELEMENT(x,l)/ELEMENT(x,2)) ELEMENT(x,4)]), x, 
[a,m,l,0]),4),m)) 

For those of you who are still here I am going to unpack now a little present as a reward: A truly 
powerful powermod- function! One day in the not too distant future the programmers of DERIVE will 
be merciful to those of us who are in desperate need of a fast and efficient built-in routine to compute 
large powers a n mod m, but until then take this one: 

PM(a,n,m) == IF(n < 0, PM(INV(a,m), -n,m), 

ELEMENT (ITERATE( [IF(MOD(ELEMENT (x,3),2) = 1, 

MOD(ELEMENT(x,l) ELEMENT(x,2), m), ELEMENT(x,l)), IF(ELEMENT(x,3) > 0, 
MOD(ELEMENT(x,2) A 2,m), ELEMENT(x,2)), FLOOR(ELEMENT(x,3),2)],x, 
[l,a,n]),l)) 

There are a lot of useful applications for this utility function. For instance, it can be used to emulate 
RSA public-key cryprosystem with realistic numbers, as it was shown in [2]. Furthermore many 
primality tests require the computation of large powers a n mod m. Let me conclude with some 
examples of this kind: 
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Some examples of primality testing using the subsequent powermod-function PM 

INV(a, m) := 

#1: If GCD(a, m) = 1 

MOD(ELEMENT (ITERATE (I F(MOD(ELEMENT(x , 1), ELEMENT (x, 2)) = 0 


2)), ELEMENT (x , 4 ), ELEMENT (x, 3) - FLOOR (ELEMENT (x , 1)/ELEMENT 

PM (a, n, m) := 

If n < 0 

#2: PM(INV(a, m) , -n, m) 

ELEMENT (ITERATE([IF(MOD(ELEMENT (x , 3), 2) = 1, MOD(ELEMENT (x 


IF(ELEMENT (x , 3) > 0, MOD(ELEMENT (x , 2) A 2, m) , ELEMENT (x, 2)), I 
1 0 ^ 1 00+9 has no prime divisor < 1000 , as it is shown by 
PRIM(n) := 

If NEXT_PRIME(n - 1) = n 
#3: n 

1 


100 


#4: 


GCD 


10 


+ 9, 


1000 

n 

n=2 


PRIM(n) 


= 1 


Even so, it is not a prime by Fermat's theorem: 

100 100 

#5: IF(PM(2 , 10 + 8, 10 + 9) = 1, 1, 0) = 0 

needs 0.031 sec - in 1994 2.0 sec. 

And with the built-in NEXT_PRIME function ... 

100 100 

#6: IF(NEXT_PRIME(10 + S) = 10 + 9, 1, 0) = 0 

0.041 sec against 3.7 sec. 

On the other hand ... 


100 100 

#7: NEXT_PRIME(10 ) - 10 = 267 

100 100 

#&: IF(PM(2 , 10 + 266, 10 + 267) =1, 1, 0) = 1 

It passes also the Miller-Rabin test for base 2: 


' 100 
10 + 266 100 

#9: PM 2, , 10 + 267 

2 j 


100 

-10 - 267 = -1 


Johann Wiesenbauer provided a DERIVE 6 file for his Titbits 2. You can download it together 
with all other files connected with contributions of DNL#13. 
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list2 := ITERATE 


#19: 


IF GCD 10 


plOOo] = 1, APPENDS 


rx p. 

. 2 _ 


+ 2 


x, [[], 1], 499 


1 


#20: list2 : 

= [9, 

37, 

39, 61, 63, 

99, 

103, 

111, 

123 , 

133, 

151, 

159, 

177 , 

207, 

211, 

223, 

229, 

253, 

261, 

267, 

271 , 

273 , 

277 , 

289 , 

291, 

321 , 

331 , 

363, 

369, 

391, 

411, 

433, 

453, 

463, 

499, 

511, 

529, 

531, 

547, 

559, 

579, 

583, 

589, 

597, 

609, 

613, 

627, 

639, 

643, 

649, 

663, 

667, 

669, 

679, 

691, 

697, 

711, 

721, 

729, 

733, 

741, 

753 , 

757 , 

763, 

771, 

781, 

783, 

793, 

811, 

817, 

831, 

849, 

853, 

867, 

889, 

919, 

921, 

937, 

939, 

949, 

957, 

973 , 

981, 

991, 

993] 









#21: DIM(1 ist2) = 85 

#22: plOOOO := TTCk , k, VECTOR (p r -im (k) , k, 10000)) 

#23: plOOOO : = 


5949 06 79 5799986 35868564 02 2242 50344 56 5968 03 2244 06 79 3 279383 09 70391614 04056 38 38 34 34 7448 3 32 2055 54 5213 05 5814077 2894 046169 7 0074 7:L864099413- 
17 55 06658216194 7914 0098 74 3288 33 51062 728565 33978719 76967 7794 73 72S7419461E9992 24 56614 2 568 26 7097 2713 3515885 581jS 78614846 298 2514 28 0965124- 
95 59 05459895 36 2 39 739285 29 05 0987 086 07 3243 5109 79404 549164 0561Q117 75 74461721E42601179 3 28 33 3049 78 7849113 799967 376074 07179 049797 301069941- 
80691956740996319139767941335979146924614249205S620511jS209408]j815447356600355657S589364425892Q173935047432096609314309876S9142265284- 
9886 0747 05 04 744 264182113137 76 07 3 36 5498 02908749 23 081028763 23 56 26416 5194 05 55 529817 22 39 2 396747971129 369 2697 59 28 0075 08 06 0848026 239665966- 
66156 3617910149 3 59 54 7 2 04 2 038 3 5 5 5 3 002138 29134 029 56 08 26 3 06 2 Q13 269 056 7414109 066 7 0029 59 2 74 5134 78 56 2 03 346 748 58 02 59 03 54 5 3988 713 38664 34 2 56 5- 
6 2999666615 074 2 2 08 5 7 5 02 56 3 011098 7 59 3488 5 0149 7 74 08 344 364 3 2469 07 5 014 3 00761413 7 74 317 599849 26 248 2 7 2 5 5 398 2 5 2 7 014 00124 5 31069 09866444612 BIB- 
7683 77 79 05 01248666415 31488351016 07 03 0024 791*75 0613 2 0795 278343 7248 2 36 57 77 50579092103 59154059 254083839 5 2408614436158429 3967 388997 25 04 7- 
764551147675224303476105002718019447014918746493066055578 515146 77316 523034726 044 3 56 368 3 3 7 78 711315 024 56 012 3 5 7 05959816 327 0169883153 218- 
5 03 24 06 5 29 296468 014888481S699916 2 2441124 59 264916 2 5 24 5 5 7156 3 74106 2 58 7 2 7 7 7 36 3 7 2 09 000712 59 0866 059 0314 0499 76 2 06 281064 7 03 28 5 7 288 79 581159 5- 
964 CCL313 3 5 2 7 244 0349 2 04 312 05 07 59 7 05 7 0046 3 3 7 2 01110088 386 7 318 588 59414 7 2 39 26 2 38 05 512 7 09896864 5 2 08 7 5 74 049 7617898 566981218 7813139 000044 06 3- 
411548841804721531710646179536615178801439888671896228240817295393922057018201449223072236273458767074652512060052626222363111614497- 
8 59418 34 002 79 58 54996108 3 2 2186 2 3 04 2 5179 217618 3 015 261517129 28 79 09 289 3 3 798 3696 78844 5 76 216 7 3 5 3 78 5 5 5 2 3 39 059 3 59 7 315 5 7 2124 76 049 3 3 7 5 3 36 3412 0- 
4 561£ 7 03 24 7 36 715 2610615 2 348 3 5 0415 56 5 51569 2 3 2 5 5 7 3 2 3 06 04 07648 3 2 5 56 5 048 7124 78 5 2 7 58 3 3159812 04846 341£ 5 7 09 5134 56 74 7 018 2 3 3 049126 2 28 517 2 7 2 7 0- 
3 7 299 211391244 34 059 29 36174 2 21176 78126 05861*816 2 3 5 008115 24 08 2134 7 0101717 3 2 04 7 599886 3 2 2 24 09 7 7 7 31044 3 02 74 2115698115 3126 54166 3 212124 24 5 7- 
161574538634380394023168771522864681988916 0363 26 06 5 78 778749 29 24 03 5 7179 268783 20819741346 37 0140264519 215 36 5763 38 24 33 22 26 74006 5108 38 05 5- 
35119025415817887 07 56 52 758045 539 56 5968 0445 521263 38 3 302314 344662 0488899 36 50859153 58 5 38 0124 24 054057 33 08417 33 04 78 0482402 03 2416 3107 23 713- 
2 28494 3 088 3 7 2 7 3 5 5 2 39 7 04116 5 56 046 7 00749 5 3 00068 5 215 7 06416 084917 5 3 3 2 7 5817 215 02 512136 3 74 7 0549 78108 049103 7 088 3 7 2 09 2 2 03 08 5 2 3 79 7 3 008861596 5- 
767962389150115866366580330193859432999862851517233780964251173790562077974559634515899719889044664493190077601924672112096921288946- 
917 04 3 5 364819813 04 09 3 3 3996 5 34 2 5 08 78 389 06415 2 0548 2898 38 2 58412 346448 7 599691248 59168 2 7 004 21988 7 03 38 3 3 599 7 2 34819 03489 316488 764 0217 009966- 
86817244 736947119 28 56 29 0493 55809 02 72 0617919362829 2 015441744 55 2168 286 5417 3568 79 24 7291984 5588 3895 7916 0017037 22 55 28421691544 28 08139396 5- 
41702893 73 2917 0629580544995 25 5496264 551916 5884 2064 247 04 74 2615 789714617 29 71001949 76 7 3 083 35 2685 05414284 088 52812561126368 5 7344 57 56 02928- 
33 38999598 0698 74589 324383 2547 0071662434761929586 017 35 33626025 5598 5817012671512 24 2044615 79 78 2815468515 04092 5292817115 3 7 76966 764617 751- 
2 07 5 09 712109 515 2 7 3846 3 78 2 5 09 2 2 217 08 015 39 3 564 3 2 09 79 3 5 7698156 26 2 039 02946 07 7 7 05 024816 2 5991944 29941464 2489 2 09 5 216115 2 024 344 007 0596849 7 02- 
7076 2248 24 389964 02 597 508914 0695 78 36740989 312 39 096 3 09126 038 099 09 016 7 21197 361416668 5644817138 07815 56117 02 583 2312 710039 0413985 38 03 53 517- 
95 26 77 32 24 07 306 089511.7612 72821915284 5768124159 089 5 74045 75 710389 3617 398 3449 289126 574141159 37469027405 74 724013 5948 2497 5 02 06 7814 596008 5- 
5 7 7 2 5 0798 3 5 2126 2124 29448 5 3 31549644108444134 3866446 38 OS 7696 7613 3 7 028179 3 088 54 04 3 02886 58 5 7 3 2813 028 76 74 2 3 2 3 3 366 5158 7699 5 3 2 24 0686 3 717 514- 
48 29 03 516154 515 32 0542 74 75 2105688 3157669581111915 2 24 71578078 2684906 7 26078 04 26 5064 578 56958124 7796 20559 3441336042 2545 024 54646980009177 2- 
9 0888 7 26 099 3 5 59 7489 268 03 7 3 3413 7121412 34 5 05986915158 0212 2913613669446 3 7 0154 2 21546 2 2 5 59 74 014 056 2 5 5 3669 38 74 7 2 3 7 003 74 068 54 2 217 34 06 215 389- 
8 516 7 5 2 5416 7 2 56 38 58 0266 2 005 5 0048 00124 2 78884 7 315 217 3694 3 28 02469 7 3 5 2 7113 21074 28 2 0469 717 21448 5108869 2 7 7 2696 5116 2 2 3 5 00964 5 2415 244968 5 004- 
326453057611380122048887987244539567207333743647215113231043534965839270947604956877850310506878523001614337579343647743516735315593- 
948553895278516786301660848157173547793336058716488376311646305506131123276703531083537913044515349040175385838806346081134261562257- 
57314 283 269948 216 017 294 09 5 0028 59 51584 2 5 7 748116 741716763753452 711113 0468 710 


0.5 see!! 


4297 

#24: plOOOO := 5.949067957-10 

fermatCv, a, b) := (ITERATE (["IF (PH (a, b + v - 1, b + v ) = 1, APPENDS , [v 1) , x ) , x + 1], x , [[] , 1] , 
#25: L x|2 x|2 1 L x|2j 1 2 J 


DIM(v))) 

1 


100 

#26: fermat(list2, 2, 10 ) = [267, 949] 


Dr. Wiesenbauer's email address, if anyone wants to contact him is 

J. Wiesenbauer 6 tuwien . ac . at 

[1] Nobauer W. - Wiesenbauer J., Zahlentheorie, Prugg Verlag, Eisenstadt 1981 

[2] Heugl H. - Kutzler B. (ed.), DERIVE in Education - Opportunities and Strategies, Chartwell- 
Bratt Ltd, 

1994, p 51 -61 
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The Girl with Long Hair 

J.M. Cardia Lopes, Porto, Portugal 


Sometimes DERIVE gives unexpected but funny pictures with PLOT. Please see what DERIVE returns when 
we try to display the direction field for the equation together with two particular solutions: 


y_iZ +2j + i = o 


DIRECTION_FIELD 


4-y 


4 'j 

2-x , x, -1, 1, 10, y, -1, 1, 10 

x > 





J-LDTS 

Plot IS 
t0= 5 

, , , 4 ■ y 1 _ , 4 

■/yl 1 - A 

yi 1=. 5 
y2'=l 
yi2= 
y3' = 

yi?= 

■j4 1 


v 2 ' <Q= 


SELECT DNE iST-DRDER FUNCTION ONLY 


deSoloe^y 1 - ~~ + 2x+-^- = 0,x,yj 

y = G1 x 4 + x 2 + 1 
soloe(y = @1 -x 4 + x 2 + 1 ? Gl) | x = 1 and y =► 

ei = -i 

soloefy = @1 -x 4 + x 2 + 1 ? Gl) | x = .5 and y^ 

G1 = -12.00 


■ *x A 4+x A 2+l 01 > I x= ■ 5 and v= . 5 

IAIN RAD AUTD DE 3/30 


4 2 

x - x + y - 1 


#3 : SOLVE 


4 2 

x - x + y - 1 


o, y 


4 2 

(y = - x + x +1) 


4 2 

#4: TABLE(- x + x + 1, x, -3, 3, 0.001) 


4-y 4 

#5 : DS0LVE1| + 2-x + , 1, x, y, 0.5, 0.5 

x x 


4 2 

12 -x -x + y - 1 


= 0 


#0 : SOLVE 


4 2 

12 -x - x + y - 1 


0, y 


4 2 

= (y = - 12-x + x +1) 


4 2 

TABLE(- 12-x + x + 1, x, -3, 3, 0.0CCL) 


f^JSlz^lT^cefRe^aphltl^hlD^JlVIlC 1 

t t t t / / " \ \ \ 

f t f y y / / -ia \ 

f f I f ! I f fTfr 

/ / / y / / y/y ypi 

\ t / -~ \ \ \ \ \ \ 

y vu \ U 

M M M 

\k \\\ \ \ \ \ \ \ 

/ 7 / / 7 7 7f 7 f 1 

y y y y y y / y / i 

i y y y y y 1 i f i 

i i / / / / A i i i 

i i i i i i li i li i 

1 I ''V l V V V 1: \ 1: 
1 1 \ \ \ \ \ \ \ \ 

1 Yi 1 A \ \ \ \ \ 1 

1 I I \\\ \ \ \ 1 1 

1 lil ill 1 1 1 1 1 

MAIN RAD AUTD DE 


# 7 : 


